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Production of p-wave pions in nucleon-nucleon collisions 
is studied according to an improved power counting that em- 
bodies the constraints of chiral symmetry. Contributions from 
the first two non- vanishing orders are calculated. We find 
reasonable convergence and agreement with data for a spin- 
triplet cross section in pp ppi^", with no free parameters. 
Agreement with existing data for a spin-singlet cross section 
in pp pmr'^ constrains a short-range operator shown re- 
cently to contribute significantly to the three-nucleon poten- 
tial. 



The use of (approximate) chiral symmetry of QCD 
to determine the form of the low-energy effective La- 
grangian has proven to be a powerful aid to the un- 
derstanding of strong interaction physics |^] . ft has 
long been known that the use of chiral symme- 

try for pion-nucleon (irN) scattering leads to a qualita- 
tive understanding of the pion-ranged part of the three- 
nucleon force, believed to produce important effects in 
nucleon-deuteron (Nd) scattering 1^ and few-nucleon 
bound states ||^. Yet, discrepancies between theory and 
experiment (for Ay in Nd scattering and for excited 
levels in bound states ||] ) remain which have been widely 
attributed to unknown three-nucleon forces. A novel 
three-nucleon force, expected on the basis of power count- 
ing arguments, involves the exchange of a pion between 
one nucleon and two others interacting via short-ranged 
forces This force can indeed affect Nd scattering at 
a currently observable level, and thus potentially resolve 
the remaining discrepancies |^]. ft depends on a pion- 
two-nucleon interaction of a form determined by chiral 
symmetry, but strength determined by parameters, di 
of Eq. (||), not fixed by symmetry. We argue here that 
the production of p-wave pions in nucleon-nucleon (NN) 
scattering offers a unique opportunity to determine di 

In the last few years, the various NN — > NNtt reac- 
tions have been studied both experimentally and theoret- 
ically |T^, with a focus on near-threshold energies. The 
first high-quality data concerned the total cross section, 
and most theoretical analyses have concentrated on 77 
0.4, a region dominated by the Ss state. (Final states 
are labeled by LI with L and I being the relative angular 
momentum of the nucleon pair and the pion with respect 
to the two-nucleon center of mass, respectively; 77 is the 
maximum pion momentum in units of the pion mass, 
m^). Many different mechanisms are expected at these 
kinematics: heavy meson exchanges |12 , (off-shell) pion 



rescattering , excitations of baryon resonances [|5| , 

and pion emission from exchanged mesons p^ ]. 

The pion dynamics are largely controlled by chiral 
symmetry constraints, and the hope that the use of Chi- 
ral Perturbation Theory (xPT) would yield insights led 
to the use of tree-level yPT to calculate the cross sections 
close to threshold |l7|-|2l[|. Ref. emphasized that the 
diverse contributions to the Ss final states can be or- 
dered in powers of \/ raT^ IMqcd, where Mqc d — 1 GeV 
is the typical QCD mass scale. The implication of this 
relatively large parameter is that loop diagrams enter at 
next-to-leading order in s-wave pion production. Thus 
a test of the convergence of the series is hindered. We 
shall show that such difficulties are not present for the 
case of TT production in p- waves (r\ ^ 1), because the pro- 
duction proceeds through leading-order operators better 
determined from other processes. 

Our arguments rely on the use of symmetries. One 
may obtain the results of QCD by using the most general 
Lagrangian involving the low-energy degrees of freedom 
(pion TT, nucleon iV, and delta isobar A) which has the 
same symmetries as QCD. These are approximate chiral 
symmetry, parity and time-reversal invariance. Chiral 
symmetry plays a crucial role in low-energy processes be- 
cause it demands that, in the chiral limit where the quark 
masses go to zero, the pion interactions contain deriva- 
tives, which are weak at small momenta, Q. Because the 
quark masses are small, any non-derivative pion interac- 
tions are also weak. Although the nucleon mass rriAr is 
not small, it plays no dynamical role at low energies. The 
delta isobar can be excited, but its mass difference to the 
nucleon, bm = rriA — rnjsr, is not large. For processes in 
which Q ~ it is convenient to introduce the "chiral 
index" of an interaction A = d + | — 2, where d is the 
number of small-scale factors, that is, derivatives, tUt^, 
and dm; and / is the number of fermion field operators. 
Chiral symmetry implies that A > Our interac- 

tion Lagrangian is given, using an appropriate choice of 
fields, by the expressions P3,M 
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where C4 = C4 + . The terms denoted by • • • include 
Hermitian conjugates, s-wave ttN scattering terms and 
terms of higher powers in pion fields. Our principle aim is 
to determine the parameters di = 0{1/ f^MgcD), which 
determine the desired three-nucleon force. Here we fix 
signs by taking gA = +1.26 in the chiral limit. The Ci 
have been determined from TriV scattering at tree level 

^^(tree) ^ g^^-l ^(tree) ^ 2.25 GeV"! ||) 

as well as to one- loop order (C3 = —5.29 GeV~^ and 
^{loop) ^ g GeV"^ !§). (Since we treat the delta 
isobar explicitly, we subtract its contribution from these 
values of q |^.) As will be established below, up to 
next-to- leading order, the di, which support only S ~* Sp 
transitions, are the only undetermined parameters in p- 
wave pion production. 

The next step is to extend the power counting of 
Ref. jl7j to the region 77 ~ 1, where the outgoing pion 
has energy lu = 0(771,^) and momentum \q\ — ©(mTr), 
and the two nucleons in the final state have momentum 
\p'\ — ©(m^) and total energy p'° = 0{m'^/mN)- The 
unique difficulty of using xPT for pion production is that 
the entire pion energy is supplied by the relatively large 
momentum of the initial nucleons, \p\ = 0{y^mppiri^). 
Note that the non-relativistic approximation holds, as 
p'^/Smlf m^/mAT <C to^t ~ p2/2mjv. 

The scales of momenta and energy are not the same, 
so it is simpler to count powers of the small scales in 
time-ordered perturbation theory. Equivalently, one first 
integrates over the time component of loop momenta 
in covariant diagrams. In this case, an intermediate 
state is associated with an energy denominator a 
loop with a (3'^/(47r)2, a spatial (time) derivative with Q 
(i?), and a virtual pion vertex with from wave 

function normalization. For N, E ^ Q'^/mpf, for A, 
E ~ Q^/ruN + Sm, and for tv, E ^ \/ Q'^ + m^- 

Final-state interactions (FSI) are those which occur 
after the emission of the real pion. In this case, the nu- 
cleons have typical Q m.^. The energies of interme- 
diate states containing a tt or A can be ii^ ~ m^r, but 
otherwise E ~ ni^/mN- The sum of "irreducible" sub- 
diagrams where all energies are 0{mTr) is by definition 
the NN potential, which is then amenable to a xPT ex- 
pansion. The sum of "reducible" sub-diagrams produces 
the final-state wave function {ipf). 

In contrast, all intermediate states occurring before the 
radiation of the real pion are characterized by loop mo- 
menta ~ ^Jrof^rn^. For these kinematics we find that any 
additional loop requires at least i) one more interaction 
— pion exchange or shorter range — with an associated 



factor no larger than 1//^; ii) a volume integral with an 
associated factor of {^^JrfC^^miqY j \ and iii) an addi- 
tional time slice. If the additional time slice cuts a pion 
line, a factor of 1/ ^m-^ratq comes in, and the overall ex- 
tra loop factor is at least 1 (V™^™ ^ — 1 = 2lh^ 
that is, a suppression by two powers of the expansion pa- 
rameter. If the additional time slice does not cut a pion 
line, a factor of appears, and there is a relative 

enhancement of i/TOjv/m^. Integrals over two-nuclcon 
states typically also have enhancements by factors of tt 
from the unitarity cut. Thus we resum those diagrams 
that differ by the addition of interactions between the 
initial nucleons (ISI) , and the effects are contained in an 
initial state wave function IV'i)- 

These considerations yield a pion production ampli- 
tude T = j-iAi). Both the kernel K and can 
be obtained from the chiral expansion, but the currently 
available jV'i) do not yield an accurate fit to the measured 
NN scattering phase shifts. Therefore we use a phe- 
nomenological coupled-channel {NN, NA, AA) model, 
CCF of Ref. H, fitted to NN scattering. 

The leading contributions to p-wave production are 
displayed in Fig. |l]. At lowest order (0(1), apart from 
overall factors) there are contributions from the direct 
production off the nucleon and off the delta, where all 
vertices are from (Fig- 0i, ii)- At next-to-leading 
non- vanishing order (OimTr/mN)) there are four types of 
contributions. First, there is a recoil correction to the di- 
rect production. Second, there are rescattering diagrams 
that proceed through the seagull vertices in C^^^ pro- 
portional to 1/4/^ (Galilean correction to the Weinberg- 
Tomozawa term), C3, and C4 (Fig. |l|iii). Third, there 
is a rescattering through the Weinberg- Tomozawa term, 
where the primary production vertex is proportional to 
the external pion momentum. Fourth, there are short- 
range ■7t{N'^N)'^ interactions proportional to di and ^2 
(Fig. I^iv). Diagram iv) and most of the rescattering 
diagrams contribute to charged-pion production only. 
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FIG. 1. Lowest-order contributions to p-wave production. 
Diagrams at C(l) are (i,ii), and of 0(m,r/'Tijv) are (iii,iv). 
A solid (dashed) line denotes a nucleon (pion), and a double 
line a A. Interactions from £^'^^ (-C'^'') are denoted by a dot 
(circled dot). Diagrams with a A in the final state are also 
included. 

With our theory in place, we consider the available 
pion production database. This has been enriched re- 
cently by very accurate determinations of spin observ- 
ables at 0.5 ^ 77 ^ 1 for ppn'^ , pp — + pmr^ [|o) , 
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and pp dn^ . It is useful to describe the total cross 
section in terms of components '^^'^^am, where S is the 
initial NN spin with projection m along the direction of 
the incoming momentum. The ^'^"^^(t,„ can be expressed 
as linear combinations of the total cross section and the 
double polarization observables Actt and AaL 

In order to test convergence for the p-wave produc- 
tion, we need an observable where the lowest contribut- 
ing partial wave is p and the initial and final nucleons are 
not both in S states. Such an observable exists, namely 
the '^CTi cross section in neutral-pion production with Pp 
as the lowest partial waves contributing. While the ra- 
tios between double polarization observables and the to- 
tal cross section, AcrT/ctot and AtTi/atot, have recently 
been accurately measured at lUCF the total cross 
section is known to a much lesser accuracy (see the com- 
pilation of Ref. 1^). To determine the error of the total 
cross section, we simply take the total spread of the data 
as the error band. We defer a more detailed analysis un- 
til it can benefit from the soon-to-be-available much 
better data. 

In p-wave production the lowest-order loop contribu- 
tion enters one order higher, at 0{{mT^ /miq)^), than the 
rescattering terms of Fig. ^ and are ignored. Besides the 
coupling constants of the pion to the baryon fields, the 
only parameter that enters is C3. We will use both values 
given above to get an estimate of loop effects on the final 
result. 

80 




FIG. 2. Chiral perturbation theory predictions for ^ai in 
the reaction pp — » ppn^. Lowest order (long dashed line), 
lowest order plus recoil contribution (dot-dashed line), and 
next-to-leading order using 03'°°^' (solid line) and ci^^""^^ (dot- 
ted line) are shown. Data are from Refs. ps] , ^ . 

The predictions of xPT are compared to the data in 
Fig. ||. Up to values of 77 ~ 0.7 the data is well described. 
Deviations at higher energies might be due to higher par- 
tial waves entering, and/or to higher-order p-wave con- 
tributions. In any case, we see that sub-leading correc- 
tions are smaller than leading contributions throughout 
the range 77 ^ 1. 

We next consider the amplitude for the ^Sq ^ {^Si 



Di)p transition, denoted oq, which has recently been ex- 
tracted from the reaction pp pmr^ [Q . The loop cor- 
rections are again expected to be small, but the number 
of rescattering diagrams is larger, since isospin-odd op- 
erators (the recoil correction to the Weinberg- Tomozawa 
term as well as the C4 term of Eq. (H)) enter. The strik- 
ing feature of oq is that interactions proportional to the 
di's also contribute. Because there seems to be reason- 
able convergence in the p waves, we assume that they 
can be reliably computed and that we can attribute any 
deviation between theory and experiment to the effects 
of the terms involving the coefficients di. The contact 
interactions enter as the linear combination di + 4c?2- 
Thus there is one unknown parameter to be fixed by the 
data. On the basis of dimensional analysis we expect 



d = Udi+M2 



with S = 0(1). 



Our result for qq is shown in Fig. |[ We find a de- 
structive interference between direct nucleon and delta 
contributions that makes oq small and more sensitive to 
sub-leading terms. For the Ci parameters we employ the 
values extracted from the tree level fit to ttN scattering 
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We use dipole form factors; to make contact 



with Ref. we employ cutoff parameters A = 1 GeV 
for diagrams containing pion exchange and A = for 
the contact interactions. The result for (5 = is not in 
disagreement with data, whereas a value of (5 = 1 leads 
to a serious disagreement with experiment. In Ref. |^ 
d — —0.2 was shown to yield an important contribution 
to Ay in Nd scattering at energies of a few MeV. Using 
d = —0.2 here is also consistent with the pion production 
data. 
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FIG. 3. ao of pp np-K in chiral perturbation theory. 
The different lines correspond to values of the parameter re- 
lated to the three- nucleon force: S = 1 (long dashed line). 
(5 = (dot-dashed line), S = —0.2 (solid line), and 5 = — 1 
(short dashed line). Data are from Ref. |34| . 

In contrast to ^cri, the result for is quite sensitive 
to the cutoff parameter used in the rescattering contri- 
bution, because the momentum range scanned by the C4 
term is quite large. For example, our results for ao can 
vary up to a factor of 2 if the corresponding cutoff pa- 
rameter is increased to 2 GeV. The cutoff sensitivity is 
not a serious difficulty because it also occurs in calcula- 
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tions of three-nucleon forces. From the viewpoint of an 
effective field theory this can be simply understood: the 
large momentum pieces of the loop integrals involved in 
the evaluation of the C4 contribution can be absorbed by 
a counterterm, namely d2- Thus, the cutoff dependence 
of C4 directly translates into a scale dependence oi d2- A 
reasonable phenomenological estimate should follow from 
using the same cutoff and parameter set in both calcu- 
lations. On the experimental side, it is clear from Fig. 
^ that a reduction of the uncertainty in the data would 
allow a stronger constraint on S. We find this a strong 
motivation to the continuation of the existing program 
on pion production. 

We have shown that there is convergence in p-wave 
pion production, and that data on this reaction can 
be used to extract information about the three-nucleon 
force. It is clear that more accurate data would be very 
useful. In particular, the parameter d could be extracted 
and the calculation of Ref. |^ repeated to predict three- 
nucleon observables. We find it very gratifying that chi- 
ral symmetry provides a direct connection between pion 
production at energies ~ 350 MeV (lUCF) and Nd scat- 
tering at energies ~ 10 MeV (Madison, TUNL). 
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Humboldt Foundation. 
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